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C ■ Coarse-grained flux density profiles in type-II superconductors with non- 

^ • parallel vortex configurations are obtained by a proposed phenonienological 
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I ^i ' least action principle. We introduce a functional C[H] which is minimized 

^~~^ • under a constraint of the kind J £ A{H,x), where A is a bounded set. 
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t^^ ■ penetration profiles H{x, t) are derived when a changing external excitation is 

^^ . applied. Faraday's law, and the principle of minimum entropy production rate 

"j^ ■ for stationary thermodynamic processes dictate the evolution of the system. 
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in crossed field measurements. 
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The magnetization curve of type-II superconductors may display physical properties 
against the expectations of equilibrium thermodynamics. In particular, the existence of 
hysteresis and non-negative magnetic moment has been routinely observed. The Critical 
State Model (CSM), which dates back to the work by C P. Bean |[I|, has been an essential 
phenomenological framework for the interpretation of the aforementioned experimental facts. 
The following prescription was given: External field variations are opposed by the maximum 
current density Jc within the material. After the changes occur Jc persists in those regions 
which have been affected by an electric field. 

Currently, the irreversible properties of superconductors are well understood in terms of 
the vortex flux line lattice (FLL) dynamics in the presence of pinning centers. Within the 
framework of self-organized extended dynamical systems one can conceive the CSM as the 
competition between a repulsive vortex-vortex interaction and attractive forces towards the 
pinning centers [0]. This results in metastable equilibrium states for which the gradient in 
the density of vortices is maximum, corresponding to the critical value for the macroscopic 
current density Jc- At the macroscopic level, one usually makes the assumption that the 
rearrangement to new equilibrium states is instantaneous whenever the system is perturbed. 
As a matter of fact, when the FLL is unpinned by an external drive, a diffusion process 
is initiated, which is characterized by a time constant Tf ~ ^oL'^/pf (pf stands for the flux 
flow resistivity and L is some typical length of the sample). Thus, the previous hypothesis 
corresponds to neglecting Xf as compared to the excitation typical period. 

A serious limitation of Bean's model is that one can just apply it to lattices of par- 
allel flux tubes. However, a wealth of experimental phenomena is related to interactions 
between twisted flux lines. Macroscopically, if vortex crossing is present, J develops full 
vectorial character and the condition \J\ = Jc (or 0) does not suffice. To the moment, 
several phenomenological theories are at hand, which allow to deal with such cases. Among 
them we want to detail the work by Clem and Perez- Gonzalez 0. These authors have 
developed a model (double critical state model or DCSM in what follows) which includes 
current density components perpendicular and parallel to the local magnetic induction B. 



Their corresponding critical values Jci_ and Jc\\ are respectively associated to depinning and 
flux-cutting phenomena. Metastable distributions of B are obtained by the critical-state 
principles | J±| < Jc±i \J\\\ < Jc\\i appropriate E{J) laws and the Maxwell equations. Finite 
element based models are also available ||^ , which allow to compute the critical state profiles 
for non-ideal geometries. 

Here, we will show that Bean's simplest concept of opposing external field variations with 
the maximum current density is still valid for multicomponent situations; the sign selection 
for one dimensional problems will become a particular case of finding the adequate direction 
of J by way of a minimum principle. Variational approaches based on the minimization of 
the free energy for calculating the magnetic properties of type-II superconductors have been 
applied before. In Ref. a numerical method is presented which allows computing the 
current distribution for the Meissner state in finite cylinders. Also, in a previous work 
we showed that a powerful generalization of variational calculus, the optimal control (OC) 
theory [0, provides a very convenient mathematical framework for critical-state problems 
in superconductors. Generally speaking, the OC tools may be fully exploited in physical 
theories that include limitations in the form of inequalities. In Ref. |0 , the magnetostatic 
energy was minimized under the restriction \J\ < J^. That principle may only be applied to 
the initial magnetization curve, as hysteretic losses cannot be accounted by a thermodynamic 
equilibrium model. Nevertheless, it was suggested that a functional, most probably related 
to changes in the magnetic field vector should allow dealing with the full problem. In this 
letter we show that the OC tools may be used to predict the irreversible quasistationary 
evolution that takes place. Eventually, the theory will be applied to the experiments of a 
superconductor in rotating magnetic fields and the magnetization collapse in crossed field 
measurements 0. 

The basic relations of coarse-grained electrodynamics in the case of type-II supercon- 
ductors read as follows. As time-dependent phenomena are involved, one must incorporate 
Faraday's law f^odtH = —V x E (we have used B = fioH, which means that reversible 
magnetization is neglected in this work), as well as an appropriate E{J) characteristic for 



the superconductor, where J is to be obtained from Ampere's law V x H = J. 

We will assume a discretization scheme in which H^ stands for the magnetic field intensity 
at the time layer n6t. This procedure permits posing the minimum principle as a tractable 
boundary value problem for ordinary differential equations. In order to gain physical insight, 
we will infer the CSM equations after considering some aspects of the more familiar eddy- 
current problem in normal metals {E = pj). The successive field profiles in a magnetic 
diffusion process may be obtained by the finite-difference expression /io(-f?n+i — H^)/5t = 
— pV X (V X i^n+i); which defines a differential equation for Hn+i- Notice that for each step 
one can identify it as the stationarity condition for the functional 

Chi[Hn+l\ = /^O / l-f^n+1 — -f^iil + St E ■ J = / jFjj+i , 

Jn Jn Jn 

where Q stands for the sample's volume and the dependence of the second term on //n+i is 
implicitly assumed. In fact, the Euler-Lagrange equations which describe the stationarity of 
Cm, i.e.: dJ-'^i+i/dH^+ii = i9^.[i9jFn_j_i/9(9^.i/n+i,i)] can be checked to produce the aforemen- 
tioned expression. We call the readers' attention that Cm should not be mistaken for the 
action in the classical theory of fields S = J JC dV dt. As an additional advantage of using 
Cm, we get a clear physical picture of the underlying series of quasistationary processes. 
Notice that Cm holds a compensation between a screening term and an entropy production 
term. In fact, under isothermal conditions one has S = E ■ J jT . Thus, a perfect conductor 
would correspond to the limit 5^0^ -f^n+i -^ -f^n- On the opposite side, non-conducting 
media would not allow the existence of screening currents (otherwise S -^ oo) and, thus, 
//n+i "will be solely determined by the external source. In the case of type-II superconduc- 
tors, the critical state arises from the flux flow characteristic, which, in the isotropic case, 
can be written as E = pf{J — J^) (or if J < Jc). Thus, the external drive variations are 
followed by diffusion towards equilibrium critical profiles in which J equals Jc- If the relax- 
ation time Tf may be neglected (or equivalently pf -^ oo) the superconductor will behave as 
a perfect conductor for J < Jc and as a non-conducting medium for J > J^. In the light 
of the previous discussion, the evolutionary critical state profiles can be either obtained by 



using Maxwell equations and a vertical E{J) law or the principle: 

In a type-II superconducting sample Q with an initial field profile H^^ and under a small 
change of the external drive, the new profile H^+i minimizes the functional 



C[#n+i(f)] = i /|#,+i-#,|2, 
^ JO. 



for the boundary conditions imposed by the external source, and the constraint V x /fn+i G 
A(i?n+i,x). 

For simplicity, we will use the isotropic hypothesis |V x iJn+i| < -^c(|-f^n+i|) hereafter, 
i.e.: A is a disk. This will provide a nice agreement of our simulations and the experimental 
facts. However, anisotropy can be easily incorporated, for instance by choosing A to be an 
ellipse or a rectangle (DCSM case) oriented over different axes. Although isotropy would not 
seem to be justified according to the underlying physical mechanisms of flux depinning {Jci_) 
and cutting (Jc|j), it may be supported by other reasons. As a matter of fact, an average 
description seems adequate for highly twisted soft FLLs for which flux cutting phenomena 
are much more effective than for rotating rigid parallel sublattices |T^ . 



In order to see how the OC machinery arises, let us consider an infinite slab of thickness 
2a in a field parallel to the faces {YZ plane) and take the origin of coordinates at the 
midplane. By virtue of the symmetry, we can restrict to the interval < x < a. Along 
this work, we will use a Kim's model type Jill dependence of the critical current density 
Jc{H) = Jco/i^ + Hj Hq), which incorporates the microstructure dependent parameters Jq^ 
and Hq. For convenience we will express x in units of a, H in units of Hq, and J in units of 
Ho/a. Then we can state Ampere's law, together with the critical current restriction in the 
following manner 

— * 

f{Hn+i,u,x) 



dx l + \H^+i\ 

Above we have introduced the dimensionless constant (3 = Jc^a/ Hq and the so-called control 
variable u, which is a vector within the unit disk D. Notice that, by construction, one has 
u J- J. Thus, we have the state equations for the state variables Hj^^i{x). 



Next, we require the minimization of the functional C[H^^i{x)] constrained by the state 
equations. Just in the manner of Ref. [^, Pontryagin's maximum principle can be used to 
solve the OC problem. On defining the associated Hamiltonian 

the optimal solution (i.e.: functions if*^^(x) and 'll*{x) minimizing C and satisfying the state 
equations) fulfils the Hamiltonian equations 



dx ^' ' dx "+^'^ "'' ^'dH*, ' 



'n+l,i 

together with the maximum principle condition 

ueD 
In the case under consideration, the control variables must take the form u* = p*/p*, and 
this leads to the system 



(la) 



-^n+l,i - -^n,i + -fjl /i ,"rr* \9 • (^^) 



dx p* 1 + i^n+i 



A part of the boundary conditions required to solve this system of differential equations 
is given by the external field values at the surface iJ*^^(l) = H[l, {n + l)6t]. The remaining 
boundary conditions will be supplied, at every instant, according to the particular situation: 
(i) If the new profile matches the old one at a point < x* < 1, i.e.: H*^^{x*) = Hj^^x*), 
these are the extra boundary conditions, x* can be determined by additional equations 
derived from the minimum cost requirement. In fact, one can prove that a free final param- 
eter X* leads to the algebraic condition 7i{x*) = 0. (ii) If the new profile holds a variation 
which reaches the center of the slab, the full arbitrariness of i/*^^(0) supplies the so-called 
transversality conditions for the momenta: p*{0) = 0. 

Notice that the physical counterpart of the result |'ii*(a;)| = 1 is |J| = Jc{H). We 
should emphasize that this condition and the distribution rule for the components of J are 

6 



determined by the selection of the control space A. For instance, | J\ is no more fixed when A 
is a rectangle. Instead, the optimality produces a vector leaning on the boundary, matching 
the evolution predicted by the DCSM. 

Eventually, the critical state profiles will be solved by integration of the set of Eqs.(l). 
Below, we apply the method to the rotating and crossed field experiments. For definiteness, 
we choose /? = 1. 

— * 

First, we consider the solutions for H^_^^ in a field cooled sample, which is subsequently 
subjected to a surface field rotation in the manner Hsit) = Hs{0, sin as, cos as), where as = 
ujt. On neglecting the equilibrium magnetization contribution, the slab holds a nonmagnetic 
initial state of constant profile (0,0,-^^5). Successive profiles of the penetration field i?n+i 
were obtained by means of Eqs.(l). Figure |I| displays the main features of the calculated 
magnetization process in our system. During the initial stages of rotation the magnitude of 
H (upper panel) is decreased towards the center of the slab in a flux consumption regime. 
Simultaneously, the angle of rotation a respect to Hs{0) follows a quasilinear penetration 
profile (lower panel). As rotation is continued the field modulus penetration curve develops 
a V-shape, which neatly defines a decoupling point xq. Thereafter, the curve essentially 
freezes and the fiux density modulus becomes stationary. On the other hand, the rotation 
angle variation is blocked in the range < x < Xq. Further changes of Hs will only affect 
a{x) for xo < a; < 1. In particular, after decoupling occurs, the external drive variations 
induce a conventional critical state behavior for the profile a{x). For instance, one can 
observe the expected effect of rotation reversal after one turn is completed (see the inset in 
the lower panel). Eventually, the outer region will be responsible for the hysteretic losses as 
the inner part contains an inert magnetic ffux density distribution. 

The phenomenological matters described above have been experimentally observed by 
Cave and LeBlanc |^ and reported by Clem and Perez-Gonzalez |^ from the theoretical 
point of view. However, we want to remark that DCSM model contains critical slopes for 
the field modulus and rotation angle, and the appearance of a decoupling point is somehow 



forced. Our variational principle allows showing that, even for the isotropic hypothesis, in 
which no such a priori condition is introduced, the optimal process itself produces the actual 
current distribution and generates decoupling. Thus, the behavior observed in Fig.|l| is more 
related to the imposed boundary conditions than to the particular region A in use. 

Next, we concentrate on the so called magnetization collapse, which can be observed in 
crossed field measurements. To our knowledge, Ref. [§] displays a remarkable manifestation 
of this effect in high Tc superconductors. We have calculated the field penetration profiles 
for a zero field cooled sample to which a constant excitation Hzs is then applied. This is 
followed by cycling stages of the other field component on the surface Hys- Fig.^ displays the 
predicted magnetization curves. As usual, we have defined M = {H{x)) — Hs- The major 
loop (sequence OABC) displays the observed experimental features: (i)My shows a nearly 
conventional GSM profile, except for the fact that the loop is not closed (see A and C). (ii) 
Mz irreversibly collapses as Hys is cycled. Both effects may be easily explained in terms of 
the predicted penetration profiles. In the insets we show the evolution corresponding to the 
branch A — i> B. For illustration, we have also included a few profiles associated to the first 
field reversal steps in the branch B — i> C. Notice that Hy follows the typical GSM pattern, 
whereas {H^) continuously increases. Physically, this behavior must be related to the most 
effective mechanism for minimizing J |//n+i — H^l'^dx. The current density component Jz 
dedicated to reduce the imposed field variation |/7n+i,y — -^n,yP is privileged. Then we have 
Jj^ ~ near the surface owing to the restriction on \J\, and this leads to the flattening of Hz- 
We have also simulated a minor loop, which corresponds to a partial penetration regime. It 
is noteworthy that quite different behaviors can be observed depending on the applied field 
amplitudes. 

In summary, we have presented a phenomenological critical state model that general- 
izes the minimal tool proposed by Bean to systems of twisted vortex configurations. Our 
theory may be used to understand the experimental features of rotating field experiments. 
Within the isotropic hypothesis, the model also provides a straightforward explanation of 
the observed magnetization collapse, for which a merely approximate justification was avail- 
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able [|T2|. Although this work has been developed for an isotropic relation of the kind 
J < Jco/(l + H/Hq), several extensions may be implemented if dictated by the physics of 
the problem. These include the effect of equilibrium magnetization by means of an appropri- 
ate B{H) relation, the selection of the model Jc{H) and the use of anisotropic control spaces. 
Another important issue would be the incorporation of time relaxation effects by means of a 
finite flux flow resistivity. This can be accomplished by using the entropy production term 
E ■ {J — Jc) in the functional. 

The authors acknowledge financial support from Spanish CICYT (project MAT99-1028) 
and from DGICYT (project DGES-PB96-0717). 
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FIGURES 
FIG. 1. Magnetic field penetration profiles for a field cooled slab under successive rotation 

steps for the surface vector Hs- Panel (a) displays the modulus consumption towards the center of 

the slab (x = 0) as well as the tendency towards a stationary frozen V-shape. Panel (b) displays 

the rotation angle a respect to the initial constant profile. Subsequent to the appearance of the 

decoupling point xq, which has been marked on both graphs, the evolution restricts to the range 

xo < X < 1. The inset shows the calculated angle profiles upon rotation reversal. H has been used 

in units of Hq, x in units of a and a is given in radians. 



FIG. 2. Evolution of the magnetization components in a simulated crossed field experiment for 
a zero field cooled superconducting slab. Subsequent to the application of a constant surface field 
Hzs, the other component Hys was cycled either in a major loop (sequence OABC) or minor loop. 
Several magnetic field profiles, corresponding to the magnetization process have been included 
in the insets. Full symbols have been used for the curves corresponding to the points A and B, 
continuous lines for a selection of intermediate profiles and dashed lines for the initial steps in the 
branch B^C. Hy is plotted within the axis range (-0.6,0.6), H^ within (0.25,0.65), and x for (0,1). 
All the quantities are in dimensionless units as defined in the text. 
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